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ABSTRACT 

We study generation of magnetic fields by the Biermann mechanism in the 
supernova explosions of first stars. The Biermann mechanism produces magnetic 
fields in the shocked region between the bubble and interstellar medium (ISM), 
even if magnetic fields are absent initially. We perform a series of two-dimensional 
magnetohydrodynamic simulations with the Biermann term and estimate the am- 
plitude and total energy of the produced magnetic fields. We find that magnetic 
fields with amplitude 10~ 14 — 1CT 17 G are generated inside the bubble, though 
the amount of magnetic fields generated depend on specific values of initial con- 
ditions. This corresponds to magnetic fields of 10 28 — 10 31 erg per each supernova 
remnant, which is strong enough to be the seed magnetic field for galactic and/or 
interstellar dynamo. 

Subject headings: magnetic fields — cosmology: interstellar medium — supernova 
remnant 



1. Introduction 

Magnetic fields are ubiquitous in the universe. In fact, observations of rotation measure 
and synchrotron radiation have revealed that magnetic fields exist in astronomical objects 
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with various scales: galaxies, clusters of galaxies, extra-cluster fields, etc (for a review, see 
e.g., Widrow 2002). The observed galactic magnetic fields have both coherent and fluctuating 
components whose strengths are comparable to each other (Fosalba et al. 2002; Han et al. 
2004). Conventionally these magnetic fields are considered to be amplified and maintained 
by dynamo mechanism. The coherent component in a galaxy is expected to be amplified by 
galactic dynamo, while the fluctuating component may be amplified by interstellar dynamo 
driven by turbulent motion of the interstellar medium (ISM) (Balsara et al. 2004). However, 
the dynamo mechanism itself cannot explain the origin of the magnetic fields: The seed 
magnetic fields are needed for the dynamo mechanism to work. 

Thus far various mechanisms have been suggested as a possible source of the seed field. 
They can be classified into two types, astrophysical and cosmological origins. Here we 
concentrate on the former (for the latter scenario, see e.g., Lesch & Chiba 1995; Davis 
& Widrow 2000; Bamba & Yokoyama 2004; Berezhiani & Dolgov 2004; Yamazaki et 
al. 2005; Takahashi et al. 2005). Basically, the astrophysical magnetogenesis invokes the 
Biermann mechanism (Biermann 1950) which is induced by the electric currents produced 
by the spatial gradient of the electron pressure not parallel to that of the density. This is 
a pure plasma effect so that there is no need to assume unknown physics as is often done 
in cosmological models. Because the Biermann mechanism requires the non-parallel spatial 
gradient of the pressure and density, some non-adiabatic process is necessary to produce 
deviation of the equation of state from the polytropic one. The strong magnetic fields in 
high-redshift galaxies (Athreya et al. 1998) imply that the significant amount of the seed 
magnetic field should be generated at early stage, e.g., epoch of cosmological reionization or 
protogalaxy formation. For instance, Gnedin et al. (2000) studied the generation of magnetic 
fields in the ionizing front, and found that magnetic fields as high as ~ 10~ 18 G in virialized 
objects can be generated. Kulsrud et al. (1997) showed that magnetic filed of ~ 10 -21 G 
can also be generated at shocks of large-scale structure formation. 

In this paper, we investigate magnetogenesis at smaller scales. Specifically, we study 
the amplitude of the magnetic field produced by the Biermann mechanism when the shock 
waves of the supernova explosions of the first stars are spreading throughout the ISM. The 
primordial supernova explosions are expected to take place effectively, since the initial mass 
function (IMF) of population III stars should be substantially top-heavy (e.g., Abel et al. 
2002). We perform a series of the two dimensional magnetohydrodynamic (MHD) numerical 
simulations in which the Biermann term is included. We also discuss whether they can be 
the origin of the cosmic magnetic fields. Consequently, we find that the spatially-averaged 
amplitude of the produced magnetic field in virialized objects reaches ~ 10~ 16 G, which is 
much greater than those expected from cosmic reionization and large-scale structure for- 
mation. Thus the supernova explosions of the first stars can be effective sources for the 
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seed magnetic fields. Although the situation considered here is somewhat similar to that of 
Miranda et al. (1998), they assumed multiple explosion scenario of structure formation and 
considered explosions of objects with mass > 1O 6 M at z > 100, which is clearly unrealistic 
in the context of the current standard model of structure formation. 



2. Numerical Simulations 
2.1. Numerical Method 

The basic equations for our numerical simulations are as follows: 
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where p, v, e, P, and B are density, velocity, total energy, pressure, and magnetic field, 
respectively. Here the total energy of the gas is defined as, 



where we fix the value of the adiabatic index 7 to 5/3, which is valid in case of the non- 
relativistic gas. The last term of the right hand side of equation (4) is the Biermann term: a 
in equation (4) is the so-called Biermann coupling constant defined by a = m p c/e(l + x) ~ 
10~ 4 G • s, where m p , e, and x being the proton mass, electric charge, and ionization fraction, 
respectively. Although the gas temperature just before the star formation is rather cool 
T ~ 200K (Abel et al. 2002; Bromm et al. 2003; Omukai & Palla 2003), UV radiation 
from the first stars ionizes the surrounding ISM (Freyer et al. 2003; Mori et al. 2004). 
Thus we set x = 5/6 assuming nne/nn ~ 0.1. The radiative cooling is represented by A in 
equation (3) and we use a cooling function derived by Raymond et al. (1976). When parcels 
cool below 10 4 K, an artificial heating rate proportional to the density, T in equation (3), is 
used. The constant heating coefficient is set so that heating balances cooling at the ambient 
density and temperature. Although the cooling function and heating rate in the primordial 
gas are not clear so far, they are not important for the adiabatic expansion phase which we 
concentrate on. 
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We solve the above equations by the two-dimensional MHD code in the cylindrical 
coordinates (r,z,(f>) assuming axial symmetry around the symmetry axis (z). The code is 
based on the modified Lax-Wendroff scheme with an artificial viscosity of von Neumann 
and Richtmyer to capture shocks. The numerical scheme was tested by comparing known 
solutions which have been obtained either analytically and numerically. Specifically, we 
checked the code by comparing with [1] the adiabatic SNR with the Sedov solution (Sedov 
1959), [2] the spherically symmetric SNR (Slavin & Cox 1992) and [3] MHD shock tube 
problems (Brio & Wu 1988). The analytical solution of [1] is reproduced within 5% relative 
error, and the oscillations behind the shock front are well suppressed. 

The artificial viscosity is added to smooth the jump at the shock front in equations (1) 
- (4). For all the quantities U = (p, pv, E + B 2 /87r, B), we added an extra term to express 
artificial diffusion and solved equations like 

^ = V(cVC/), (6) 



where we take 
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with A = Az = Ar and A, = 4. 



In all the computations, grid spacings are chosen Ar = Az = 0.1 pc. For example, 
numerical domain covers a region of 130 pc x 130 pc with 1300 (r) x 1300 (z) mesh points 
in our fiducial model (see §2.2). Even if the resolution is doubled, the total energy of magnetic 
field is not changed 4 times at t — 10 5 yrs. We begin the simulation by adding thermal energy 
of E = 10 53 or 10 52 erg within the sphere of 2pc in radius. 



2.2. Results of Numerical Simulations 

As the bubble expands, the ejected gas interacts with ISM and shock wave is formed. 
In the shocked region the gas is heated non-adiabatically, which is a necessary condition 
for the Biermann mechanism to work. Here the structure of the interstellar environment is 
important because it affects the density and pressure profiles of the shocked region which is 
directly related to the Biermann term. We assume inhomogeneous ISM with average density 
^ism = 0.2cm~ 3 . This is roughly consistent with the situation discussed in Bromm et al. 
(2003). The scale length of the density and the amplitude of the inhomogeneity are poorly 
understood now and we assume inhomogeneity with the scale length A = lpc and density 
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variation 0.2 x 2 ±1 cm~ 3 , which are similar values to those in our galaxy. Within the variation, 
the amplitude of density is given at random and the distribution is smoothed numerically 
to perturb. This is our fiducial model for the ISM. As we will show in the next section, the 
amplitude of the produced magnetic field is sensitive to scale length A, while the average 
density and the amplitude of the density variation is rather unimportant. Thus we consider 
several different models in addition to the fiducial model: Specifically, we vary the mean 
density (1 x 2 ±1 cm~ 3 and 10 x 2 ±1 cm -3 ) and the scale length (3pc and lOpc). 

As for the explosion energy of the supernova, we adopt -Esn = 10 53 erg for the fiducial 
mode. This explosion energy corresponds to stars with mass 250M© which explodes as a 
pair-instability supernova (Fryer et al. 2001). Also we consider a model with = 10 52 erg 
as a variation. 

Figure 1 shows the contours of the gas density and amplitude of the produced magnetic 
field at the end of the adiabatic expansion phase t = 1.26 x 10 5 yrs for the fiducial model. 
The radius of the bubble is about 125 pc and turbulent motion is induced in the shocked 
region due to the inhomogeneity of the ISM. The amplitude of the magnetic field is about 
10~ 14 G for the central region and about 10~ 17 G just behind the shock. The total magnetic 
energy inside the bubble is about 10 30 erg. 

We have also checked the case of homogeneous medium to test robustness of our compu- 
tation. We find that SNR generates the magnetic fields with amplitude ~ 10 -19 G behind the 
shock front in average, and it is ~ 10~ 3 times smaller than that in inhomogeneous medium. 
Therefore, in case of homogeneous medium, we estimate a numerical error of the amplitude 
of magnetic field as ~ 10~ 19 G. It should be noted that the produced magnetic field has only 
toroidal component because axial symmetry was assumed. 

In Figure 2, we show the time evolution of the total magnetic energy for various models. 
The behaviors are qualitatively similar for all the models. The robust knees around 10 3 years 
in Figure 2 come from a formation of adiabatic shock front: It corresponds to when SNR 
shifts from free expansion phase to Sedov phase. The total magnetic energy at the end of the 
adiabatic expansion phase is larger for models with smaller scale length of the ISM density. 
This tendency will be confirmed by an order-of-magnitude estimate in the next section. For 
models with large average ISM density or small explosion energy, the total magnetic energy 
is smaller because the bubble is smaller than the other models. Although there are many 
uncertainties in initial conditions, the generation of magnetic field with the total energy 
10 28 — 10 31 erg appears to be robust. 
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3. Analytic Estimates of Magnetic Fields in Primordial Supernova Remnants 

and Implications for the Seed Magnetic Field 

To understand the result of numerical simulations, in this section we perform an order- 
of-magnitude estimation of the strength of the magnetic field produced by the Biermann 
mechanism. 

The amplitude of the magnetic field produced by the Biermann mechanism can be 
estimated from the Biermann term in equation (4) 

Vp x VP , . 

^Biermann ~ Oi At, (8) 

P 

where At ~ 10 3 yr is a characteristic timescale the Biermann mechanism works. Taking 
the characteristic pressure to be the ram pressure of the gas, P ~ P ram = p ^ ubble , and the 
characteristic velocity of the bubble v bubble ~ 10~ 3 pcyr -1 , we obtain 
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where L is a scale length of the pressure component perpendicular to the density gradi- 
ent. Then the magnetic energy produced by the Biermann mechanism per each primordial 
supernova remnant can be estimated as 
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which is consistent with the value obtained from our numerical simulations. 

The dependence of the total magnetic energy on several parameters can also be under- 
stood from equation (10). It is found that E B oc A~ 2 directly from equation (10). To examine 
the dependence on the other parameters, we simply assume the Sedov- Taylor solution: 

Pbubblc OC t 2/5 (P S N/niSM) 1/5 , ^bubble OC t~ 3/5 (PsN/niSM) 1/5 • (H) 

Putting these into equation (10) yields E B oc (Psn/^ism) 7 ^ 5 - Our numerical results shown 
in Figure 2 are quite consistent with these simple estimations. 

Now we estimate the spatially-averaged energy density of the magnetic fields produced 
by the first stars and consider if they can be a source of the seed fields. As for the primordial 
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star formation rate, we extrapolate the one by Pello et al. (2004) and Ricotti et al. (2004); 
p± ~ 10~ 2 M Q yr _1 Mpc -3 . Denoting the magnetic energy produced by the Biermann mech- 
anism as €sn ~ 10 30 erg, the magnetic energy density produced during the formation period 
of the first-star (r ~ 1 Gyr) can be obtained as 



M SN - 

40 (h-\( —J^ ^ ( J^-Y ( -J™} ( -L-\ er e cm- 3 ( 



where Msn is the typical mass scale of first stars that end up in pair instability supernovae, 
and / 77 is the mass fraction of such first stars: We adopt / 77 = 0.06 that was derived 
under the assumption that very massive black holes produced from first stars end up in 
supermassive black holes in galactic centers (Schneider et al. 2002). We note that the value 
is the comoving density averaged in the universe: We can convert the value to physical 
density in virialized objects (i.e., protogalaxies) as 

_ 34 / e B \ fl + z\ A ( A \ _ 3 
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where A is the density contrast. This corresponds to the mean magnetic field of B ~ 10~ 16 
G in protogalaxies, which is much stronger than expected in ionizing fronts, B ~ 10~ 18 G 
(Gnedin et al. 2000). 



4. Summary and Discussions 

We have studied the generation of magnetic fields in primordial supernova remnants. 
We have performed the two dimensional MHD simulations with the Biermann term which 
can produce magnetic field through the non-adiabatic interaction between the bubble and 
ISM, even if there is no magnetic field at first. We have found that the ISM around the 
primordial supernovae is an effective site to produce magnetic fields. The the total energy 
of the magnetic fields is 10 28 — 10 31 erg, depending on parameters adopted. Based on the 
results, we have estimated spatially-averaged energy density of the magnetic fields produced 
by the first stars during the formation period of the first stars. The averaged energy density 
is about 10~ 40 erg cm~ 3 , which corresponds to B ~ 10~ 16 G in protogalaxies at z ~ 10. This 
is much greater than expected from cosmic reionization and large-scale structure formation. 
Thus primordial supernova remnants would be a promising source for the seed fields for 
galactic and/or interstellar dynamo. 

We note that our results are based on two-dimensional MHD simulation. This mecha- 
nism may behave very differently in more realistic three dimensional case, primarily because 
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the vorticity cascade may be different. We are preparing for the three-dimensional MHD 
simulations to study this possible difference, which will be presented elsewhere. 

Although the coherence length of the seed field computed here is much smaller than 
the galactic scale, it can be amplified by the galactic dynamo to produce the coherent 
component if the coherence length is about 100 pc (Poezd et al. 1993; Beck et al. 1994; 
Ferriere & Schmitt 2000), which is a typical size of supernova remnants. It might be also 
amplified by interstellar dynamo to produce the fluctuating component (Balsara et al. 2004). 
While it is beyond our scope to discuss the relation between the produced magnetic field 
and galactic/interstellar dynamo, we plan to investigate the evolutions of the seed magnetic 
fields computed here as a result of the dynamo processes in the large scale of galaxies and 
cluster of galaxies. This will be presented in the forthcoming paper (Hanayama et al. in 
preparation). 

As observational signature, a proposal to detect seed magnetic fields was given by e.g., 
Plaga (1995). If there exists intergalactic magnetic fields produced by the primordial SNR, 
an arrival time of high energy gamma-ray photons from extragalactic would be delayed by 
the action of intergalactic magnetic fields on electron cascades (Ando 2004). The magnetic 
field as weak as ~ 10~ 24 G would be detectable if the delay time comes into reasonable range, 
a few days. Therefore the mechanism of magnetic field generation proposed in this paper 
might be tested in the future high-energy gamma-ray experiments such as GLAST. 

We would like to thank Kohji Tomisaka for invaluable comments and advice on writing 
the manuscript. We are grateful to Kanako Sugimoto and Hiromi Mizusawa for useful 
comments. Numerical computations were carried out on the VPP5000 supercomputer at the 
Astronomical Data Analysis Center of the National Astronomical Observatory, Japan. We 
also thank Ryoji Matsumoto and Takaaki Yokoyama for a contribution to the calculation 
code, CANS (Coordinated Astronomical Numerical Software). K. T., K. K., M. O. and K. I. 
are supported by Grant-in-Aid for JSPS Fellows. 
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Fig. 1. — Contours of the gas density (left) and the magnetic field (right) at the end of the 
adiabatic expansion phase (t = 1.26 x 10 5 yr) for the fiducial model. The amplitude of the 
magnetic field is about 10 _14 G at the central region and about 10 _17 G just behind the shock. 
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Fig. 2. — Time evolution of the total magnetic energy produced by the Biermann mechanism 
for various models. 
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